Roughening Transition in a One-Dimensional Growth Process 
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A class of nonequilibrium models with short-range interactions and sequential updates is presented. 
The models describe one dimensional growth processes which display a roughening transition be- 
tween a smooth and a rough phase. This transition is accompanied by spontaneous symmetry 
breaking, which is described by an order parameter whose dynamics is non-conserving. Some as- 
pects of models in this class are related to directed percolation in 1-1-1 dimensions, although unlike 
directed percolation the models have no absorbing states. Scaling relations are derived and compared 
with Monte Carlo simulations. 
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The morphology of growing interfaces has attracted 
much interest in recent years jl). Many growth pro- 
cesses of two-dimensional surfaces exhibit a roughening 
transition, from a smooth phase with finite width to a 
rough one with diverging width. A question of interest 
is whether a one dimensional (Id) growing interface with 
short-range interactions and unbounded noise, can ex- 
hibit a roughening transition It is well known that 
in thermal equilibrium no such phase transition can take 
place as Id interfaces are always rough. Growth processes 
far from equilibrium are, however, less restrictive and 
the question of whether they are capable of exhibiting a 
roughening transition in Id is more subtle. Most growth 
processes, such as those described by the KPZ equation 
are always rough in Id. A class of Id models, which 
have a maximal velocity by which the uppermost point 
of the surface can propagate, has been shown to display 
a roughening transition 1^,^. The existence of a maxi- 
mal velocity in these models is due to the use of parallel 
updates, and the smooth phase disappears if sequential 
(continuous time) updates are used. Sequential updates 
are a more adequate description of systems where differ- 
ent particles are not exactly synchronized. The question 
of whether a sequential update growth process is capable 
of exhibiting a transition from a smooth to a rough phase 
is still open. 

A related and more general question is that of spon- 
taneous symmetry breaking (SSB) and long range order 
in Id systems |^. Recently a nonequilibrium Id model 
with short-range interactions and unbounded noise which 
exhibits SSB in the thermodynamic limit was presented 
The model belongs to a class of driven diffusive 
systems, in which charges of two kinds are injected at 
both ends of a Id lattice and are biased to move in op- 
posite directions. The microscopic rules are symmetric 
under space and charge inversion. However, this sym- 
metry is broken in the steady state of the system where 
the currents of the two charge species are different. In 
this model SSB is a result of the conserved dynamics of 
the order parameter in the bulk (charges are not created 



or annihilated, except at the boundaries), and the exis- 
tence of open boundaries (two endpoints) at which the 
dynamics is different from that of the rest of the system 
Q. These two features create favorable conditions for 
SSB. The conserved dynamics slows down the evolution 
of the system; moreover flips between one broken sym- 
metry phase to another can originate only at the two 
boundary points, where the order parameter is not con- 
served. Simple attempts to modify the model such that 
either one or both of these features are eliminated results 
in symmetric steady states with no SSB. It would there- 
fore be of interest to examine the possibility of SSB in Id 
systems under more general conditions, namely in homo- 
geneous systems with periodic boundary conditions and 
order parameters with non conserving dynamics. 

Finally, phase transitions in homogeneous nonequilib- 
rium Id systems have usually been observed in the past 
in systems which have absorbing states (a set of states 
from which the system can not escape). The canonical 
example is the 'dry' state below the percolation thresh- 
old in directed percolation models Thus it is of 
interest to find Id models with no absorbing states that 
display a phase transition. 

In this Letter we present a class of nonequilibrium 
models with short-range interactions and sequential up- 
dates, which addresses the three questions posed above: 
The models describe Id growth processes which display a 
roughening transition between smooth and rough phases. 
This transition is accompanied by SSB associated with a 
non-conserved order parameter in a homogeneous system 
with periodic boundary conditions. The models supply a 
robust local mechanism for eliminating islands of minor- 
ity phases generated by fluctuations in the bulk of the 
majority phase. We derive some of the scaling properties 
of a particular model in this class which can be related 
to directed percolation pO|-p^. However, unlike directed 
percolation, the model has no absorbing states (to be dis- 
cussed below). The scaling predictions are compared to 
Monte-Carlo simulation results. 
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A. Model Description: The class of models is most 
simply introduced in the language of interface growth 
psf , in which both adsorption and desorption processes 
take place. In the present models, desorption may take 
place only at the edge of a plateau. For concreteness, 
we study two particular models in this class, (a) a re- 
stricted solid on solid (RSOS) version that may also be 
considered in a charged particle representation, and (b) 
an unrestricted model that may be related to directed 
percolation. Both models are on a Id lattice with pe- 
riodic boundary conditions and are defined as follows: 
Let hi be the (integer) height of the interface at site i, 
i = 1 . . . N. The interface evolves by choosing a site i at 
random and carrying out one of the two following pro- 
cesses: (a) adsorption of an atom 

hi hi + I with probability q (1) 

and (b) a desorption of atoms from the edge of a step 

hi — > min(/ii, with probability (1 — '?)/2 (2) 
hi — > min{hi, with probability (1 — q)/2 (3) 

In the RSOS version, the restriction 

\hi - hi+i\ < 1 



(4) 



is imposed at all sites. The RSOS version may be viewed 
as a driven diffusion model of two oppositely charged 
types of particles. The charges 



Ci^i+i = hi+i - hi e {-1,0, +1} 



(5) 



are bond variables and represent a change of height be- 
tween adjacent interface sites (see Fig.l). In this repre- 
sentation, the dynamical rules (|l|)-(§) correspond to ran- 
domly selecting two neighboring bonds and performing 
the following processes with probabilities as indicated on 
the arrows: 
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In both RSOS and unrestricted models, when q is 
small, a smooth phase is maintained. In this phase the 
interface displays a small concentration of short-lived is- 
lands, and its average velocity, v, is zero in the thermo- 
dynamic limit. As q increases, adsorption increases and 
typical islands grow, until, above a critical value gc, is- 
lands merge and full new layers are completed, giving the 
interface a finite growth velocity. Thus, when q is small, 
a local mechanism that eliminates islands is present in 
the model : an island is formed with boundaries that are 
biased to move towards each other (due to desorption 
from the island edges). The evolution of a large island 
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FIG. 1. (a) Typical configuration of the interface, is 
the fraction of sites at height k above the minimal height in 
the configuration (here no = ni = |, n2 = |). The average 
island size grown on top of level k is It- (b) Mapping of the 
configuration of Fig. la to the charged-particle representation, 
along with a site coloring, as described in Sec. |^ 

for (7 < is illustrated in Fig. 2. It is seen that the is- 
land shrinks, ensuring the stability of the smooth phase. 
This behavior is typical of islands of all sizes, except the 
very largest (i.e. a complete layer). The rule that no 
holes can be formed in a completed layer (00 -/-> — h 
in the RSOS version), preventing it from dividing into 
shrinking islands, is essential for obtaining the smooth 
phase. To demonstrate the existence of the roughen- 
ing transition, we carried out Monte Carlo simulations 
of both models. In this Letter we present some of the 
results obtained in this study. A more detailed account 
will be published elsewhere ||lj]. The phase transition 
takes place at q^ = 0.189 ± 0.002 for the RSOS model 
and = 0.233 ± 0.001 for the unrestricted model. The 
interface width is defined by the standard deviation of the 
height distribution w = [N-^ ^-(/i, - N'^ J^i h-ifV^^- 
We find that starting from a flat interface, w rises as 
w ~ t^^^ for short times, saturating for large t at w ~ 
where N is the lattice size. At q > qc, the numerical re- 
sults are consistent with the KPZ exponents § X = | 
and z = |, indicating a rough interface. Below qc, w sat- 
urates at a finite value independent of N , showing that 
the phase is smooth. The critical behavior at q — q^ 
shown in Fig. 3, is 



w - log(7V) 



(7) 



In the following we discuss the symmetry breaking which 
takes place for q < qc- We also discuss the relation to 
directed percolation and the critical behavior near qc- 

B. Spontaneous Symmetry Breaking: To demonstrate 
some of the model's properties, it is convenient to con- 
sider the RSOS model in the charged particle representa- 
tion (Eq. 1^, Fig. 1). At (7 < (7c, the charges are arranged 
as closely bound H — dipoles. At q > q^ the dipoles be- 
come unbound, and the fluctuations in the total charge, 
measured over a distance of order N, diverge with TV. 
Thus the transition is manifested in correlations between 
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FIG. 2. Monte-Carlo simulation of the RSOS model for a 
system of size 600 at g = 0.130 where = 0.189±0.002. Each 
configuration is a row of pixels, with sites at even and odd 
heights represented by black and white pixels, respectively. 
Configurations at intervals of 7 moves per site (sweeps) are 
shown (time advances downwards) up to 2100 sweeps. At an 
early time, a large island of size 400 is introduced. The is- 
land shrinks and disappears, illustrating the mechanism that 
insures long-ranged order at g < Qc. Note that small islands 
and islands within islands are continually generated by fluc- 
tuations, and are washed away. 

charged particles rather than in their density. The sym- 
metry breaking which takes place in this model is best 
seen by introducing a coloring scheme by which each of 
the sites between the charged particles is colored in one 
of two colors, such that the two sites adjacent to a -I- or — 
particle have different colors, and the two sites adjacent 
to a particle have the same color (Fig. lb). Every move 
in the dynamics corresponds to a local rearrangement of 
the charged particles and site coloring. Under the model 
dynamics, any configuration of charges and coloring can 
evolve to any other allowed configuration of charges and 
coloring. Thus the model has no absorbing states. How- 
ever, at g < gc, a typical configuration displays unequal 
concentrations of the two colors. As the system evolves 
the configurations flip between majority colors in a time 
scale which was found to grow exponentially with the sys- 
tem size |l4|. Thus, though the dynamical rules are sym- 
metric with respect to the site colors, the system spon- 
taneously selects one of two colors as a majority color, 
breaking the symmetry. The system in the phase space 
of charge configurations and colorings is ergodic at any fi- 
nite size, but becomes non-ergodic in the thermodynamic 
limit when q < qc- To quantify this symmetry breaking, 
we define a magnetization-like order parameter (valid for 
both the RSOS and unrestricted models) 

1 ^ 

^ N E(-l)'' (8) 

which can be envisaged by considering the two colors as 
'up' and 'down' spins. The order parameter is clearly not 
conserved by the dynamical rules. In the smooth phase 
{q < Qc)i (M) 7^ in the thermodynamic limit. On the 
other hand, in the rough phase (M) — 0. Monte-Carlo 



simulations (Fig. 3) show that near the phase-transition, 
for both models, 

(lAfD-e", e = 0.55 ± 0.05 . (9) 

where e = qc — q. 

C. Relation to directed percolation: Some features of 
the unrestricted model may be related to a directed per- 
colation (DP) model ||lO|-[l^, allowing a derivation of 
several scaling properties. The occupation of the low- 
est exposed level (no in Fig. la) corresponds to the wet 
or percolating sites. The non percolating or dry sites are 
the sites where higher levels are occupied. A wet region 
may become dry both at its edges and bulk sites, while a 
dry region may shrink only at the edges. This defines a 
contact process [[l5[ , which is a sequential update version 
of a DP model ||lj]. The percolating phase corresponds 
to the smooth phase in the model obtained at g < qc. 
Thus, the occupation of the lowest level should vanish 
at the transition with the exponent (3 characterizing the 
critical behavior of the DP wet phase 

no-^e'"«, xo=/3«0.28 (10) 

The front velocity for q > qc may be related to the life- 
time of typical wet islands below the percolation thresh- 
old. These islands have a lifetime which diverges at the 
percolation threshold, with the critical exponent i^n of 
the DP coherence time |l^ . The time to complete a new 
layer is the time it takes for its missing regions (perco- 
lation wet sites) to be covered by adsorption (dry up). 
Thus the velocity v is proportional to the inverse of the 
island lifetime: 

vr^{-e)y, y = i.||«1.73 (11) 

These exponents are in good agreement with the values 
measured by Monte-Carlo simulations, xq = 0.29 ± 0.03 
and y = 1.7±0.1, as shown in Fig. 3. For the RSOS model, 
there is no direct mapping to DP, since layers grown on 
top of a dry island affect its evolution. Monte-Carlo sim- 
ulations and diagonalization of the time-evolution opera- 
tor of smaU systems {N < 12) [^,0, suggest that both 
RSOS and unrestricted models have the same exponents, 
and therefore belong to the same universality class. 

We now present a simple scaling argument for the be- 
havior of rifc, the density of sites at height k above the 
lowest exposed level. Consider first the lowest exposed 
level fc = 0. According to the DP picture, there are two 
length scales in the problem: the average size of the dry 
islands, Zq, which diverges at the transition as Iq ^ 
]l2t , and the transverse correlation length ^± ^ e^'^^ 
with i/± « 1.10 ||l^. The two length scales are related, 
for a system of size N, by the finite size scaling rela- 
tion ^0 ~ e-^fiNe"^), where / is a function satisfying 
/(s) ~ sl^/"^ for s ^ 0. Similarly, no ~ el^giNe"^) with 
g{s) s"^/"^-!- for s — > 0. At criticality, one therefore has 
no ~ Iq^ ~ N-f^/"^. 
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level is crucial for obtaining a smooth phase, and thus a 
roughening transition. The present models can be viewed 
as a hierarchy of DP-like processes, with the fc-th echelon 
DP process confined to the dry sites of the k — 1-echelon 
DP process. The critical exponents y and Xk are found 
to be related to DP exponents. It would be interesting 
to find out whether other exponents, such as 9 (Eq. H), 
may also be related to the DP problem. The models are 
easily generalized to higher dimensions, where the map- 
ping to directed percolation and the scaling arguments 
are expected to apply. It would be of interest to con- 
struct a coarse-grained field-theory |^,^ that describes 
the present class of models. 
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FIG. 3. Monte-Carlo simulation results for the unre- 
stricted growth model. The critical behavior of a quantity 
F ^ (q — Qc)" is estimated by finite size scaling, measuring 
F in systems of size A'' at the critical point q — q^, and us- 
ing F{N) ~ N""^"-^ , where v±_ « 1.10 is the critical exponent 
associated with the divergence of the correlation length (see 
Sec. Systems of sizes A'' = 2"\ m = 3 ... 10 were studied 
aX q = qc — 0.233 with 2^^ « 8 • 10^ moves per site, (a) 
Order parameter M. (b) Density of exposed sites on the low- 
est exposed level, no, and the next level n\. The bold lines 
have the slopes expected from the scaling arguments (Eqs. 
[lo| , p^ ). (c) Interface velocity v. The bold line has the slope 
expected from the scaling arguments (Eq. |l^). (d) Interface 
width, w, suggesting w ~ log(A'') (note that this graph is 
semi-logarithmic) . 

We now consider the level k — 1. One may view islands of 
sites at heights A: > 2 as growing on top of the dry islands 
of the fc = 1 level, whose typical size is Iq. Applying 
the same scaling relations and assuming that the system 
size may be replaced by Zoi we find ni ~ l^^ ~ 1^ ^^'^^ ^ 
where li is the mean size of islands of sites with height 
k > 2. Repeating this reasoning for the next levels, one 
has Uk ~ l^^ '■'^ Ik-i'^^ ■ therefore obtains 

nfe^e"^ xk = 13(13/ (12) 

Numerical simulations (see Fig. 3) yield xi — 0.07 ± 
0.02, as compared with xi « 0.07 obtained form Eq. 
[T^ . More extensive numerical simulations are needed to 
demonstrate the validity of Eq. ^ for this exponent and 
the exponents associated with higher levels. 

D. Discussion A class of models that describe Id 
growth processes with a roughening transition are pre- 
sented. These models provide an example of spontaneous 
symmetry breaking that takes place with a non-conserved 
order parameter in a ring geometry. The zero rate of des- 
orption of an atom who's two neighbors are on the same 
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